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Evolution of stellar structure in the Small Magellanic Cloud 
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ABSTRACT 

The projected distribution of stars in the Small Magellanic Cloud (SMC) from the Magellanic 
Clouds Photometric Survey is analysed. Stars of different ages are selected via criteria based 
on V magnitude and V — I colour, and the degree of 'grouping' as a function of age is stud- 
ied. We quantify the degree of structure using the two-point correlation function and a method 
based on the Minimum Spanning Tree and find that the overall structure of the SMC is evolv- 
ing from a high degree of sub-structure at young ages (~10Myr) to a smooth radial density 
profile. This transition is gradual and at ~75 Myr the distribution is statistically indistinguish- 
able from the background SMC distribution. This time-scale corresponds to approximately 
the dynamical crossing time of stars in the SMC. The spatial positions of the star clusters in 
the SMC show a similar evolution of spatial distribution with age. Our analysis suggests that 
stars form with a high degree of (fractal) sub-structure, probably imprinted by the turbulent 
nature of the gas from which they form, which is erased by random motions in the galactic 
potential on a time-scale of a galactic crossing time. 
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1 INTRODUCTION 



The majority of stars form in clu stered environment s (e.g . 
IClarke. Bonnell & Hillenbrand l200(j| and lLada & Ladi] 120031 
henceforth LL03), but only a small fraction (< 10 per cent) of 
stars ends up in bound star clusters which survive the embedded 
phase. It was noted much earlier that almost all sta rs in the Galac- 
tic disc form in (unb ound) OB associations (e.g. Roberts] [l 9571 : 
iMiller & Scald fl979l) and that residual gas expulsion is the most 
proba ble explanation ("e.g. lHillslll980l : Lada. Margulis & De arborn! 
1 1984b . LL03 introduced the term 'infant mortality' for the rapid 
dissolution of embedded clusters at young (few Myrs) ages. 

More recently, the infant mortality scenario has also been 
used to explain the strong drop in the age distribution, dN/dt, 
around ~ 10-20 Myr of clusters in the Antennae galaxies a nd M5 1 
dFall. Chandar & Whitmore1l2005l and iBastian et al .1120051 respec- 
tively). IPellerin et alj d2007h introduced a new approach for study- 
ing infant mortality, by comparing the spatial distribution of stars 
of different stellar types in the nearby spiral galaxy NGC 1313. 
They found that O-type stars are more strongly grouped than B- 
type stars, which is what is expected if all stars form in clusters and 
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the majority disperses on time-scales comparable to the life-time of 
O-stars (few Myrs). 

The term infant mortality is now used in widely varying con- 
texts. Firs t, the time-scales involved range from 3 Myr (LL03) up 
to a Gyr dFall et al J [20051 ; IChandar et al.ll2006l) . Second, the term 
'star cluster' is used for irregular, dynamically unmixed, groups of 
embedded stars with total masses of 10 — 100 Mq (e.g. LL03), for 
compact (few pc) and ma ssive (few times 10 5 — 10 6 Mq) clusters in 
interacting galaxies (e.g. Schweizer et al. 1996) and for large scale 
(> 100 pc) stellar associations dPellerin et alj|2007h . 

The Small Magellanic Clou d (SMC) cluster syst e m ha s been 
heavily under debate recently. iRafelski & Zaritskvl d2005l) pre- 
sented a new catalogue of SMC star clusters and showed the exis- 
tence of a declining dN/ dt. This was interpreted as infa nt mortality 
working on a 3 Gyr time-scale by 
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Chan dar et ail §006). However, 



2007) showed that the decline 



is in fact due to detection incomple teness at higher ages, later con- 
firmed by an independent analysis dde Griis & Goodwin 2008). 

This study aims at tying the above issues together to give a 
more coherent picture of the evolution of s tructures in the SMC. 
Inspired by the work of Pellerin et al. (2007), the catalogue of five 
million st ars detected in the Ma gellanic Clouds Photometric Survey 
(MCPS) (Zaritsk vetai]|2002h is used to select stars of different 
ages. In section [2] the data and the selection criteria of different 
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regions are described. In section[3]two techniques are employed in 
order to search for stellar groups and quantify their evolution with 
age. The implications and conclusion are presented in section|4] 



2 DATA AND SELECTION CRITERIA 

The MCPS covers an area of roughly 4.° 5 x 4° of the Small Mag- 
ellanic Cloud (SMC). We use the Johnson V and the Gunn I pho- 
tometric results of roughly five million stars to construct the colour 
magnitude diagram (CMD) of V vs. V — I. The photometry was 
corre cted for a Galactic foreground exti nction of Ay = 0.12 
mag dSchlegel. Finkbeiner & Davislll998l) . The resulting CMD is 
shown in Fig.Q] where the greyscale represents the number of stars 
with that colour and magnitude (dark representing higher num- 
ber). Overplotted ar e evolutionary isochrones of the Padova mod- 
els for Z = 0.004 teertelli et al.1l 19941; iGirardi et alJI 1996L l2OO0h 
w hich were converted to the UBVRIJHK photometric system 
bv lGirardi et al.l |2002h. We adopted a distance mo dulus of 18.87 
mag (~60kpc) drlarries. Hilditch & Howarthll2003l) . 

Based on the CMD we select seven non-overlapping boxes 
around the main sequence and one box around the turn-off of a one 
Gyr population, each containing 1500 stars, in a circle with a diam- 
eter of 2.7 degrees around the mean RA and DEC of the stars in the 
catalogue. This avoids the outer regions of the galaxy which have a 
higher fraction of contaminating sources. We use equal numbers of 
stars to avoid possible number dependancies in the analysis (sec- 
tion |3j- The box sizes are determined as follows: once criteria for 
the V — I colour and peak V magnitude are chosen, the stars satis- 
fying these constraints are ordered in V brightness and the brightest 
1500 are selected. The magnitude of the faintest star included in the 
final 1500 stars sample defines the lower boundary of the box. The 
boxes around the main sequence always sample stars of multiple 
ages, with the age spread getting larger for boxes lower in the main 
sequence. This means that even though the mean age increases go- 
ing down in the main sequence, we always have a contribution of 
the youngest population. Box 8 samples only turn-off stars with 
ages of roughly a Gyr. 

To quantify the contribution of the different ages to the boxes, 
we perform the following test: we assume a constant rate of star 
formation over the last Gyr (older ages do not contribute to the se- 
lected boxes on the main sequence, see Fig. []} and assign masses 
according to a Salpeter mass function. We assign colours and mag- 
nitudes to each star, based on their age and mass as derived from 
the Padova models and the adopted distance modulus. For each age 
we find the number of stars with colour and magnitude satisfying 
our criteria, allow i ng us to deri ve the mean age of stars in each box. 
See Bastian et al. U200ll2008l) for details. 

We do not take into account i nternal extinction variations in 
the SMC. lHarris & Zaritskvl J2004l) find a trend of younger stellar 
population (< 10 Myr) being more extincted than older stellar pop- 
ulations (> 1 Gyr). Especially for box 1 this could result in target 
stars being outside the box due to extinction, but we do not expect 
many sources to be shifted into our selection boxes. For the most 
extincted field in the sample of Harr is & Zaritskvl fc004h . the Ay 
distribution of the young population peaks around 0.5 mag, corre- 
sponding to E(V — I) — 0.2 mag. For most of their fields, however, 
the extinction for the young population is much lower. 

In TableQjwe summarise the selection criteria of the eight age 
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Figure 1. Colour-m agnitude diagram for stars from the catalogue of 
IZaritskv et all <2002h . The shading represents the number of stars present 
(shown as the square-root to increase the contrast), with dark representing 
a higher number of stars. The selected boxes are indicated. The numbers on 
the right-hand side of the panel refer to the base-ten logarithm of the age of 
the isochrones shown. 



boxes and the resulting minimum, maximum and mean ages of the 
stars in each box. In figure[2]we show with points the distribution 
of the 1500 stars in each age box. The lines show the Minimum 
Spanning Trees (MSTs) derived from the distribution of the stars 
(section l3.2t . 

In section[3]we will use the two-point correlation functio n and 
the statistical Q-parameter dCartwright & Whitworthl |2004|) that 
can be derived from the MST, in order to search for structure evolu- 
tion in the SMC. We compare the distribution of stars (and star clus- 
ters) of different ages to a distribution without sub-structure, which 
we will refer to as the reference distribution. This distribution has a 
smooth power-law density profile, p(r), with r the distance to the 
centre of the SMC and p(r) cx r -0,75 . Such a profile, projected in 
2D, is a reasonable approximation of the observed distribution of 
stars with V < 20. 



3 STRUCTURE EVOLUTION IN THE SMC 
3.1 The two-point correlation function 

A commonly used and relatively straightforward technique to mea- 
sure sub-structure is the two-point correlation function (TPCF), 
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Figure 2. In each panel we show the RA and DEC (relative to the mean) of the 1500 stars in each age box as dots. The mean age of the stars in the boxes 
(section|2j is indicated in each panel. The Minimum Spanning Trees (MSTs, see section |3~2l , of each sample is over-plotted with lines. 
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Figure 3. The two-point correlation function (TPCF) for stars in the eight 
age boxes. 



which has been used extensively a s a statistical m easure of large- 
scale cosmological structure (e.g. IPeeblesI 1980). We define the 
TPCF,£(r),as 



(i) 



where m is the number of stars found at a distance r from star i, 
and nl c{ is the same, but for the references distribution. For the sum 
over the reference distribution we take the average of 100 realisa- 
tions. This is similar to the procedure implemented bv lGomez et alj 
1 19931) who studied the spatial distribution of young stellar objects 
in the Taurus-Auriga molecular cloud. 

The results for all boxes are shown in Fig. [3] The slope gets 
shallower with increasing age (box number) and approaches that of 
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Figure 4. Evolution of the three structure indicators discussed in section [5] 
as a function of mean age of the eight boxes. In the top and middle pan- 
els the results of the power-law fits (equation |5J to the TPCFs are shown. 
In the bottom panel the results for Q of the clusters(stars) are shown as 
open(filled) circles. The dashed lines indicate the value corresponding to 
the reference distribution. The grey areas denote the values where the dif- 
ference between the value for the youngest age box and the reference value 
is less than 10 per cent. 



the reference distribution, that is. 
functions of the form 



1 + £(r) = 1. We fit power-law 
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Table 1. Overview of the selection criteria of the eight boxes shown in Fig.[T] 



V — I V Age [Myr] 

Box Minimum Maximum Minimum Maximum Minimum Maximum Mean 



1 -0.800 


-0.100 


8.00 


15.23 


4 


25 


11 


2 -0.500 


-0.100 


15.50 


15.93 


4 


40 


21 


3 -0.475 


-0.100 


16.40 


16.58 


4 


79 


10 


4 -0.400 


-0.125 


17.25 


17.33 


4 


126 


52 


5 -0.300 


-0.100 


17.75 


17.79 


4 


200 


76 


6 -0.300 


-0.100 


18.35 


18.38 


4 


200 


99 


7 -0.250 


-0.075 


18.90 


18.92 


4 


316 


186 


8 0.000 


0.250 


19.75 


19.76 


794 


1000 


895 



1+ ^ = a {t^) (2) 

and follow the evolution of the parameters A and B from equa- 
tion l|2} with age for the eight age boxes. This is shown in the top 
and middle panels of Fig. [4] The values evolve very quickly towards 
the reference values (A = 1 and B — 0). For the stars in box 5 
(with a mean age of ~75 Myr) the difference between the initial 
value and the reference value has already reduced by 90 pre cent 
for both A and B. This simple analysis indicates that sub-structure 
in the SMC is erased in ~75 Myr. 

3.2 The Minimum Spanning Tree and the Q-parameter 

In order to quantitatively study the distribution of groups of 
sources objective algorithms must be developed. Additionally, large 
datasets, such as the one used in the present work, require high lev- 
els of automation. One such tool which has been developed and 
successfully employed in the study of young star -forming region s 
is the Minimum Spanning Tree (MST) algorithm (Battinelli 1991). 
All points, that is, the positions of the sources, in the distribution 
are connected to their nearest neighbour, in order to form a sin- 
gle 'tree', such that the total length of all the connecting segments 
(edges) is minimised and no closed loops are formed. This is a 
mathematically unique way of describing a dataset. In Fig. [2] we 
show the MSTs for the stars in the eight age boxes. 

The MST can be used to study the distribution of 
stars a nd quantify the amount of sub -stru cture, as was done 
by e.g . ICartwright & Whitworthl d2004l) and ISchmeia & Klesserj 
(2006) who applied the MST method to study structure in (em- 
bedded) Galactic clusters. 

The Q-paramete r, which was introduced by 

ICartwright & Whitworthl J2004I) . is a measure of the amount 
of sub-structure within a spatial distribution. In short, the Q- 
parameter is defined as Q = m/s, with fh the normalised mean 
edge length of the MST and s the normalised mean correla- 
tion length. The normalisations are applied to make the result 
inde pendent of the size and number o f sources in the distribution. 

Cartwright & Whitworthl J2004I) showed that Q is a unique 
way of distinguishing between smooth radial density gradients and 
(fractal) sub-clustering. For a random 3D distributions of stars, pro- 
jected in 2D, the value of Q is 0.79. For distributions with more 
sub-structure Q < 0.79, and for smooth density profiles which are 
centrally concentrated, Q > 0.79 (see their Fig. 5). We estimate 
Q for the reference distribution through interpolation of the val- 



ues given in Table 1 in ICartwright & Whitworthl d2004l) and find 
Q rcf ~ 0.82. 

We determined Q for the stellar distribution in the eight age 
boxes and find an increase of Q with age (see filled circles in the 
bottom panel of Fig. [4). We also determine Q for the star clus- 
ters in four age bins containing equal number of cluste rs. The clus- 
ter pos itions and ages are taken from the catalogue of Hunt er et al.l 
(2003). The result is shown as open circles. The Q results for the 
stars and the clusters agree very well in the age range where they 
overlap. The ages of the clusters are known with higher accuracy 
than the ages of the stars, for which we have stars with a large 
age spread contributing to boxes 1 to 7. Therefore, the agreement 
between the Q values shows that the ages assigned by our selec- 
tion method (section O are likely to be representative. The offset 
to slightly higher Q-values for the stars is probably due to the fact 
that the stars in each box, contrary to the clusters, have always stars 
of the youngest age contributing. It indicates that the distribution of 
the stars older than the mean age affect the value of Q more than the 
distribution of the you nger stars. This is c onfirmed by simulations 
of stellar distributions (Bast ian et alj200a) . 

The arrows in the bottom panel of Fig. [4] indicate the re- 
sults for a (projected) fractal m odel with dimension 2.1 and 2.4 
( ICartwright & Whitworthll2004l) . We find that the cluster and stel- 
lar distributions in the young age boxes have a lot of sub-structure, 
with box 1 showing the same Q value as a fractal with dimen- 
sion ~2.4. This probably refl ects the structure of the tu rbulent gas 
from which they are formed (Elmegreen & Elmegreen 2001). For 
increasing age box, Q increases rapidly, levelling off around box 5 
(~75 Myr), where the value is within 10-20 pre cent that of the 
reference distribution. This is in agreement with the time-scale in- 
dicated by the TPCF It suggests that sub-structure originating from 
the star formation process is erased within ~75 Myr. 



3.3 Sensitivity of the methods 

Both the TPCF and the Q-parameter indicate that for stars in boxes 
5 and higher the distribution is essentially indistinguishable from 
the reference distribution. From a visual inspection of Fig. [2] we 
see that the distribution of stars is still evolving going from box 5 to 
box 8, while according to the TPCF and the Q-parameter these dis- 
tributions are statistically comparable (Fig.0. This could indicate 
that both our methods are insensitive to structure that is present only 
on a low signal-to-noise level. However, we note that the main body 
of the SMC also affects our results. In box 5 we can clearly recog- 
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nise the shape of the SMC, while this is less evident in box 8. This 
structure, which is of galactic scale, will increase both the value of 
Q and the slope of the TPCF. So, if the average distribution of stars 
is somewhat more centrally concentrated than our reference distri- 
bution, this could compensate the effect small structures have, such 
that the values are similar to those of the reference distribution. We 
emphasise that the absolute values of both methods are sensitive 
to edge effects and to some degree suffer from degeneracy when 
the distribution is a mixture of a pure fractal and a centrally con- 
centrated distribution. This leads us to conclude that the variations 
on the 10-20 per cent level should not be over-interpreted. How- 
ever, the general evolution of, and the agreement between, both the 
TPCF structure indicators and the Q-parameter is a strong indica- 
tion that stellar structure is erased on a short (< 100 Myr) time- 
scale in the SMC. 



4 CONCLUSIONS 

We have studied the evolution of sub-structure in the distributions 
of stars of different ages in the Small Magellanic Cloud (SMC). Us- 
ing two different methods, namely the two-point correlation func- 
tion (TPCF), and the statistical Q-parameter, we find that the stellar 
distribution evolves to the reference distribution, that is, a smooth 
power-law density profile, on a time-scale of ~75 Myr. 

Our results qualitatively agree with what lPellerin et al] f2007) 
found for stars in NGC 1313. However, we have also quantita- 
tively determined the time-scale in which 'stellar groupi ngs' dis- 
appear (~75 Myr). In the scenario of |Pellerinetai] (2007) all stars 
form in clusters/groups, which dissolve due to gas expulsion. The 
super-virial velocities of stars after gas removal are of the order 
of a few kms" 1 . Combined with the radius of the SMC galaxy 
(Rsmc — 2 kpc), the time it would take for stars that travel from 
a dissolving cluster to cover one galacti c radius is ~1 Gyr, which 
is much longer than the 75 Myr we find. IE vans & How arth (2008) 
measured a velocity dispersion, a, of 30 km s _1 for (massive) field 
stars in the SMC, which combined with the radius of the SMC re- 
sults in a dynamical crossing time, t CI = Rsmc I a — 75 Myr. This 
is the same as the time-scale for structure dispersal that we found 
in our study, making random motions in the galactic potential a 
plausible e xplanation for the erasure of structure. 

iGieles & Bastianl (2008) showed that only a small fraction 
(~3 per cent) of stars that are formed in the SMC end up in clusters 
that survive past the embedded phase, i.e. those that can be identi- 
fied as genuine star clusters in optical studies. Whether this is due 
to a 97 per cent infant mortality rate of embedded clusters or due 
to a large fraction of stars forming in the field, can not be judged 
from the current data. We do know that already at very young ages 
(~ 10 Myr) the majority of stars are not in dense clusters and these 
'field stars' still show the fractal imprints of the natal gas, which is 
being erased dynamically by random motions in the galaxy poten- 
tial. 

In an accompanying paper feastian et alj2008l) we further de- 
velop and explain the techniques and conclusions reported here, 
and compare these results to those found for the Large Magellanic 
Cloud. 
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